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Abstract 

The superconducting properties of a layered system are analyzed for the cases of zero- and non- 
zero angular momentum of the pairs. The effective thermodynamic potential for the quasi-2D 
XY-model for the gradients of the phase of the order parameter is derived from the microscopic 
superconducting Hamiltonian. The dependence of the superconducting critical temperature T c 
on doping, or carrier density, is studied at different values of coupling and inter-layer hopping. 
It is shown that the critical temperature T c of the layered system can be lower than the critical 
temperature of the two-dimensional Berezinskii-Kosterlitz-Thouless transition Tbkt at some values 
of the model parameters, contrary to the case when the parameters of the XY-model do not depend 
on the microscopic Hamiltonian parameters. 
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I. INTRODUCTION 



Theoretical description of the doping dependence of the superconducting properties of 
the high-temperature superconductors (HTSCs) is one of the most difficult problem of the 
modern condensed matter physics. Generally speaking, the complicated crystal structure of 
these materials, low-dimensional (quasi-2D) transport properties, the superconducting order 
parameter anisotropy, strong correlations and other properties result in the fact that many 
years after the discovery the microscopic mechanism of HTSC is not understood yet. 

During the last years many models which take into account some of the cuprate properties 
have been proposed. The doping dependence of the superconducting properties at T = 
in the s-wave pairing channel was studied for the 3D case in Q, Q, y] and, particularly 
for the quasi-2D case [4]. For the 2D case this problem was studied at T = ( when a 
long-range superconducting order is still possible in a 2D system (|) for the case of local 
attraction in ^1 0, 0], and for the phonon-mediated model [z| (for over- review see js}], for 
example). The ci-wave pairing along with the s-wave one at T = for the case of the 
extended Hubbard model with the next nearest neighbor attraction was studied in 
and for a 2D continuum model with short-ranger attraction and electron correlations - in 
the paper jllj]. Th properties of a model with doping dependent correlation length were 



studied recently in 



121. 



The 2D s-wave pairing at finite temperatures, when the Berezinskii-Kosterlitz-Thouless 
(BKT) transition can take place in superconducting system, was considered in for the 

case of the model with local attraction and in |l5|, |l6] for the case of the electron-phonon 
pairing. The problem of the s-wave superconductivity with fluctuating order parameter 
phase in the 3D case was analyzed in [hi \v\ . The effective action for slowly fluctuating 



c?-wave superconducting order parameter for the 2D case was also analyzed in |l8, 



19 



20 



2l|. 



However, it is known that the long-range order is impossible in the ID and 2D systems 



with an order parameter which has a continuous symmetry p. Therefore, to get real phase 
transition with a long-range order and homogeneous order parameter one needs to take into 
account the inter-layer coupling t z . The layered superconductivity is much more complicated 
since the possibility of the inter-layer fluxon and intra-layer vortex phase transitions with 
corresponding critical temperatures Tf and T v must be analyzed. It was already shown 
[22! I23I I24 . that there is only one phase transition in such a system with the critical 
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temperature T c and T v < T c < Tf ~ 8T V . The critical temperature T c is equal to T v , or, what 
is equivalent, to the temperature Tbkt of the 2D BKT phase transition at t z = 0. Then, 
this temperature value is increasing to the value Tf with the inter-layer hopping t z growth. 
In the papers (3, 23, 24, 0] the phase order parameter effective Hamiltonian was studied 
in the presence of an external magnetic field and this model was mapped on the quasi- 
2D XY-model. The XY-model parameters J\\ and J± were considered as phenomenological 
constants. It was shown that Tbkt = jJ\\ and Tf ~ 8Tbkt- 

In this paper we derive the effective XY-Hamiltonian from the initial Hamiltonian for the 
layered system of attracting fermions. In this case the parameters J\\ and J± depend on the 
bare parameters - charge carrier density, coupling, pair angular momentum, temperature and 
the inter-layer hopping. As it will be shown below, this leads to the non-trivial dependence 
of the superconducting critical temperature T c on the model parameters. In particular, in 
general this temperature is different from the critical temperature of the 2D BKT transition 
and T c < Tbkt at some values of the model parameters, contrary to the results for the case 
when parameters Jn and J± don't depend on the parameters of the microscopic Hamiltonian, 



and when the relation T c < Tbkt always holds at J± > 



II. THE MODEL AND THE THERMODYNAMIC POTENTIAL 

The model Hamiltonian for a layered superconducting system can be written as 



~\ [ dVic/Vs^tr, r 2 )^(r, r 1 )\/(r 1 , r 2 )^(r, r x )ijj j(J {T, r 2 )(l) 

where ipjcr(r, r) is a fermi-field in the with mass m and spin a =|, j, r is an imaginary time 
and j, r are layer number and intra-layer coordinate, correspondingly; t JU2 = t z (5j 2 j 1+ i + 
&J2J1-1) corresponds to the nearest inter-plane hopping. The free fermion dispersion relation 
in the momentum space has the following form 

k 2 

f (k, k z ) = h 2t z - 2t z cos(ak z ) - //, (2) 

2m 

where k is a 2D wave vector with a bandwidth W, and k z is the momentum in the inter-layer 
(z) direction, it changes in the interval [0,27r/a], where a is the inter-layer spacing; \i is the 



chemical potential of the system. In Eq.([T|) interaction V(ti,T2) describes a non-retarded 
in-plane fermion attraction. 

The partition function of the system is 



Dip^Di/je' 



with the action 



S = f dr V / d 2 r^ j<7 (T,r)d T ^ a (r,r) + H(r) 



(3) 



(4) 



To study the superconducting properties of the system with an arbitrary pairing symme- 
try the Hubbard-Stratonovich transformation with bilocal fields 0j(r, r 1; r 2 ) and <P]{t, Ti, r 2 ) 
can be applied 



cxp 



^] T (r,r 2 )^] i (r,r 1 )F(ri,r 2 )V' ii (r,ri)^ iT (r,r2) = / D(f)W(j)ex.p 



dr / d 2 r\d 



r 2 



lT,ri,r 2 J 



Let us introduce the Nambu spinor 

^■ T (r,r 



<f>}(r, ri, r 2 )^(r, r^^^r, r 2 ) - ^] T (r, rjty^r, r 2 )0j(r, r 1; r 2 ; 



(5) 



^•(r, r) 



,*](r,r) 



In this case the partition function can be written as 

Z ----- / D^DipD<fJD<f>e- S ^>^' 



(6) 



where 



Jo n J J V{n 



i\T, ri ,r 2 J 
r 2 ) 



-<9 T - r 2 



2m 



2** - m) 



%(r, r 2 )<J(r a -r 2 ) 



V I'^ 7 • r i)Tz*j 2 (T, r 2 )«J(n - r 2 ) - 5 ilj2 4(r, n, r^j^r, r^^T, r 2 ) 

- 5 iij2*n( r > r iK^ji( r > r 2)0ji(^, ri, r 2 )}, 



(7) 



where r± = |(t x ± r y ) are the Pauli matrices. 

In order to study the fluctuations of the order parameter phase and to map the corre- 
sponding superconducting effective action on the quasi-2D XY model, it is convenient to 
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make decomposition of ip^j (r, r) ^^ (r, r) on their modulus X<T t j{ T i r ) an d phase 6j(r, r), 
which as it will be shown below is proportional to the order parameter phase: 

^•(r,r) = XffJ (r,r)e^^/ 2 , 

In this case the Nambu operators are 

v& J (r ) r) = e^^/ 2 T J (r,r), 

^(r,r) = T](r,r)e-^^/ 2 , (8) 
where Tj(r, r) and TUr, r) are "neutral" Nambu spinor operators: 

Ti(r,r) = ( *f ( p r) ] ,T](r,r) = ( X ] T (r, r), Xji (r, r)) . 
VXil(^ r )y 

The order parameter can be expressed as 

&(r, n,r 2 ) = A(r, r 1; r 2 )e^( T ' ri ^ 

0t(r, rij r 2 ) = A(r, ri ,r 2 )e^^^, 

where we assume that the modulus of the order parameter A(r, r 1; r 2 ) does not depend on 
the layer index. It is also natural to assume that 

^■(r.n.ra)- A(r,r)e^ R \ (9) 

where r = ri — r? and R = (ri + r 2 )/2 are the relative and the center of mass coordinates, 

n n 

correspondingly |2lL |27| . The relation Q means that the dynamics of the Cooper pairs 
is described by the order parameter modulus the symmetry of which depends, generally 
speaking, on the relative pair coordinate and the motion of the superconducting condensate 
is described by the order parameter phase, which changes slowly with the distance and can 
be described by center of mass coordinate. In this case it is easy to obtain 

4>](t, r l5 r 2 )*t( r , n)r_^-(T, r 2 ) + *](r, ri)r + %(r, ^(r, r 1; r 2 ) 

-A^^Tt^rOr.T^r,^) (10) 
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Substituting (JHJ), © and (fTUj) into the expression for the partition function (jHJ) it is easy to 
get 

Z = J ADAD9e~^ A ' e \ 

where the thermodynamic potential is 

rP r NA(t r) 2 
/3Q(A,e)= / dr / d 2 r -TrlnG~\ 
Jo J V(r) 

N is number of the layers. The Nambu spinor Green function G can be expressed as 

g- 1 = g- 1 - e, 

where Q~ x is a part of the inverse Green's function which does not depend on the order 
parameter phase: 

^5a(n,r 2 ,ri,r 2 ) = (n, r u 3\\G~ x \r 2 , r 2 , j 2 ) 
= 5 jlh 5(r 1 -r 2 )5( y T 1 -T 2 ) 



-d Tl - t z 



V 2 

ri + 2t - (JL 



2m 



_ ^2,ii±i^( r i - r 2)5(7"i - r 2 )r z t z + 5 ili2 r x .A(ri - r 2 , r x - r 2 ). 

The self-energy E is the sum of the parts which come from the in-plane and inter-plane order 
parameter phase phase interaction E" and E^, respectively: 

E = E" + E\ 



where 



and 



Sj 1 ^(ri,r 2 ,r 1 ,r 2 ) = (n, r 1; ji|E ll |r 2 , r 2 , j 2 ) 



S jlj2 5(r 1 -t 2 )5(ti-t 2 ) 



IT 1 

^d n e h {T 1 ,r 1 )-—vi 1 e jl (T 1 ,T 1 ) 



S i"ii2( r i' r 2' r i' r 2) = <Ti,ri, ji|E- L |r 2 ,r 2 , j 2 ) 
= -^2,ii±i^( r i - r 2)5(n -T 2 )T z t z (l-exp[-iT z (9j 1 (T 1 ,ri) - 9 h (r 2 ,r 2 ).)}) 
The potential term of the thermodynamic potential is 

(3Q pot (A) =JdrJ d*r N ^* )2 - TrlnG' 1 . 
and the kinetic term can be expanded in powers of the self-energy E: 

oo 1 



(11) 



n=l 



III. THE BKT TRANSITION IN THE 2D CASE 



Let us begin with the case when there is no inter-plane coupling: t z = 0. In this case 
the behavior in each plane is independent and the system undergoes the BKT transition. 
Let us assume that the order parameter phase fluctuations are small. In this case to get 
the thermodynamic potential up to the second order in V# we neglect all the terms in (|20jh 
except n — 1, 2. Also we neglect the time dependence of 9 and the second derivative V 2 #. 
The effective potential in this case has the following structure (see, for example |8[): 

fl(A, 6) = Q pot (A) + i f d 2 r(V9f, (12) 

where 

_ f dPkdk, ( n f (k) 1 1 k 2 \ 

11 J (2vr) 3 y 4m 16m 2 T C osh 2 ( ^(k) 2 + A(k) 2 /2T) J ' 1 ' 

and the momentum distribution function n/(k) is 



„ /(k) = ! _ tanh | Vggpg (14) 
! 2T y ^(k) 2 + A(k) 2 v ; 

Free fermion spectrum £(k) in (|13j) and (|14j) is defined by (J2J) at t z = in this case. 

The minimization of the effective potential (|T2*|) at V# = with respect to the supercon- 
ducting order parameter A(k) leads to the standard gap equation: 



A(p)=/ff 1 AW tanhfvgg±ggV (Pik) , 
1 ^ 7 (2vr)3 2 A /£ (k) 2 + A(k) 2 I 2T / lP ' ; 



(15) 



(2vr)3 2v /e(k) 2 + A(k) 2 V 

The minimization of the effective potential at W = with respect to the chemical 
potential 5VL pot /5n = —vrif (v is the volume of the system) gives the equation which connects 
fi and the particle density nj in the system, or the 2D Fermi energy ep = 7mf/m: 

Ww" /(k) ' (16) 

where the momentum distribution function n/(k) is defined in (|14j) . 

To search the solutions with different angular momenta / of the pairs, we assume that 
the interaction potential has the following form: 

V(p, k) = V cos(/y? P ) cos(Zy?k)- (17) 
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Below we shall use dimensionless coupling parameter G = mV/ (2ir) for the numerical cal- 
culations. 

In the case of the interaction (|17|) the gap depends only on the momentum direction: 

A(p) = A/ cos(^p), 

where Aj is the amplitude of the superconducting gap in the case of the pair angular mo- 
mentum equal to /. The solution of the gap equation together with the number equation at 
A; = give the critical temperature of the mean-field superconducting transition Ta = T C MF 
on the charge carrier density rif. The solution of the equation 

T=|j||(A ll /*,r) (18) 

together with the gap equation and the number equation give the dependence of the crit- 
ical temperature of the BKT-transition on the charge carrier density rif. Equation (fTHjl is 
obtained by mapping (|12|) on the corresponding thermodynamic potential of the 2D spin 
XY-model. 

As it follows from the system (fTH|) and (fT£J) . the solution for the Ta and Tbkt do 
not depend on / when I ^ for the case of the simple interaction potential (|17j) . This follows 
from the fact that the /-dependence of the integral is only as cos 2 {lip) and from the identity 




where F [cos 2 (/v?)] is an arbitrary function without singularities, and I is an arbitrary non- 
zero integer number. Therefore it is necessary to analyze the solutions with I = and 

1 = 1. 

The phase diagram of the system in the 2D case is presented in the Fig.l. The temperature 
Ta is much higher in the s-channel. However, Tbkt — ep/8 in both channels at small carrier 
density. This result can be easily obtained analytically from (JT3j) and (fTHj) . 

The doping dependence of the T BKT in the cases of I = and I ^ is presented in Figs. 

2 and 3, correspondingly. The relation Tbkt — cf/8 holds up to higher values of the carrier 
density in the s-channel at fixed value of coupling. It means that the local pairs are bounded 
tighter in this case. 
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FIG. 1: Phase diagram of the 2D system in different pairing channels for the coupling parameter 
G = 1. The solid lines are Xa (the upper curve) and Tbkt for the s-wave pairing channel. The 
dashed lines are the corresponding curves for the case / 7^ 0. Here and below all quantities are 
normalized on the 2D free electron bandwidth W . 

IV. TRANSITION IN THE CASE OF COUPLED LAYERS 

Let us consider a system of coupled layers. The self-energy, proportional to the inter-layer 
coupling can be written as 

T l ± = t z T z 11i + t z Yi2 , 



where 



-<W±i<H r i - r 2)5(ri - r 2 ) cos(^ 1 (n,r 1 ) - %(r 2 ,r 2 )), 



J 2 hh 



(ri,r 2 ,ri,r 2 ) = (n, ri, |t 2 , r 2 , j 2 ) 



= <W±i<K r i - r 2)^(ri - r 2 ) sin(%(n, n) - % 2 (r 2 , r 2 )). 

Similarly to the 2D case, we assume that the phase of the order parameter changes slowly 
in the inter-layer direction. Therefore, the thermodynamic potential can be calculated up 
to the second order in (9j — 9j±i): 



Qi n = t z TTr(QT z Ei) + ^TTr{Q^Q^). 



(19) 
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FIG. 2: The doping dependence of Tbkt at / 7^ and different coupling parameters: G = 0.5 
(dash-dotted line), G = 1.0 (dotted line) and G = 2.0 (dashed line). The solid line is the function 
Tbkt = ep/8. 

The terms proportional to Sl'S -1 and S^E^ are zero due to reflection symmetry in in- 
direction. 

To map the system on the quasi-2D XY-model with the nearest neighbor interaction we 
need to obtain 

Jl 
2 



n 



kin 



J d 2 r(Vej) 2 + J z 5^(1 - cos(% - 



(20) 



This dependence comes from the first term in The second term in (flH|l is proportional 
to sin(^j — 9j±i) sin(8j — 6j±i), what is equivalent to the XY-model with the next nearest 
neighbor and next next nearest neighbor interactions. Therefore we neglect this term since 
it is of a higher order (~ t 2 z ) on the inter-layer hopping with respect to the first term (which 
is ~ t z ). However, if the coupling t z is not small this term can lead for important physical 
consequences (see for example an analysis for the 2D case Thus, the parameter J z is 



J 7 tv. 



'-rif(k) cos(ak z ). 



(21) 



Now we have obtained the kinetic part of the thermodynamic potential Qkin in the case 
of slowly fluctuating phase of the order parameter. This function is given by (120(1 . where the 
parameters Jn and J z are given by ()13j) and (|21|). Similarly to (|21|). an additional integration 
over k z must be performed in (fT3|). 
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FIG. 3: The same as in the previous Figure for the case 1 = 0. The lines for G 
G = 2.0 practically coincide with Tbkt = ejr/8. 



1.0 and for 



The effective action (fTUj) was studied in 22j,|23|, 2^1 1^ in the case when the parameters Jm 



and J z where considered independent on the fermion Hamiltonian parameters. It was shown 
that there is only one phase transition in such a system at T c which is bigger than the 
temperature of the BKT transition in the case of non-coupled layers T BKT = {ji/2)J\\. In 
the case of small coupling T c ~ Tbkt and when t z is increasing to the inter-plane hopping 
value, T c is approaching to the value T BKT = 4irJ\\ ~ 8T BKT = Tf of the fluxon transition, 
when the inter-layer order starts to take place. 

More precisely, the following expression for the effective free energy was considered 

/ dVcos^-(r) -^(r) - ^ j'" A z (v,z')dz']-E c J2 4( r )» ( 22 ) 
£o J w J(j-i)d 77 



where A(r, z) is the vector potential, 0o = hc/2e is the flux quantum, E c is the loss of the 
condensation energy in a volume ^d, £ is the in-plane correlation length, do is the thickness 
of each layer, and d(> do) is the inter-layer distance. The field Sj(r) describes vorticity of 
the lattice, Sj(r) = 1 if the vortex is present at the point, and Sj(r) = 0, otherwise. The 
length scale A e is connected with the London in-pla ne p enetration length A^ as A e = \\/ do- 
lt was shown by a renormalization group study 



24, 



251 ] that in a physical case A e 3> d the 
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self-consistent equation which describes the dependence of the critical temperature T c on 
the free energy parameters ()22j) has the form 

^ r[E c + (r/8)ln(T c /J z )] 
C ~ E c + Tln(T c /J z ) ' [Z6) 

where r = 0Q/47re 2 is connected with the BKT transition temperature as r = KTbkt- 

The comparison of the expressions ()22|) and (J2"U|) . gives the next self-consistent equation 
for the critical temperature T c , which follows from 



E c + (irJ\\/2)ln(T c /J z ) , s 

T ~47r7n— — - 117 1 v 1 ' 24) 

C ~ 11 ^ c + (47rJ||)Zn(T c /J z ) ' 1 J 

where the in-plane correlation length £ is absorbed in the parameter J z (i.e. t z {a/^) 2 — > t 2 ). 
The parameter E c actually should be renormalized by including the influence of the inter- 
layer coupling on the vortex system [2^] . It is considered here as a model parameter, which 
should be found experimentally, in particular its doping dependence should be taken into 
account. For calculation we use the value E c = 0.01W (for estimation of E c based on an 
amplitude dependent Ginzburg-Landau theory, see for example 

It is interesting to note, that in the limit of very small carrier densities, when J\\ ~ ep — > 0, 
the analytical solution for T c can be obtained T c ~ 4irJ\\ ~ ep. This is different from the 
one layer case when T c = Tbkt — £f/8, independently on the pair angular momentum I. 
However, the region of extremely low carrier densities is not interesting from physical point 
of view. 

To find the critical temperature T c one needs to solve the system of equations (fT5|). (|TB)l 
and fl2U) with functions J\\(/i, T, A(T)) and J z (n,T,A(T)) defined in (fl3jl and (l2"Tl). The 
numerical solutions show that T c < Tbkt at small carrier densities in the case of large 
values of the inter- layer hopping t z and not very strong coupling G (Fig. 4). It means that 
the dependence of the parameters J\\ and J± on coupling, carrier density and temperature 
leads to the non-trivial relation between T c and the 2D critical temperature Tbkt at some 
values of model parameters, different from T c > T B kt, as it was predicted for the case of 



fixed J\\ and J± jZ 



Qy,Q 



In general, T c is growing with the inter-layer coupling t z (Figs. 4, 5). However, in the case 
of small carrier density the critical temperature is decreasing with t z growth when / ^ 
(Fig. 4, in the / = case this effect takes place at smaller coupling G). It can be explained 
as a consequence of the fact, that the density of states on the Fermi level p(ep) at small 
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FIG. 4: The doping dependence of T c of the layered system in the case I ^ at different values of 
the inter-layer hopping and G = 1.0, E c = 0.01. The solid line is the corresponding 2D temperature 
Tbkt- The insert is the inter-layer hopping dependence of T c at G = 1 and ep = 0.05. 
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FIG. 5: The same as in Fig. 4 for the case I = 0. 



carrier densities is decreasing when system tends to become three dimensional with t z growth 
(p(ep) ~ v /ep in the 3D case and p(ep) = const in the 2D case). On the other hand, the role 
of the term ~ t\ must be studied in addition at rather large values of t z , when inter-layer 
hopping becomes of order of the intra-layer hopping, i.e. t z ~ 0.1W 7 (see again 28^). 
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V. CONCLUSIONS 



To summarize, the doping dependence of the superconducting critical temperature of 
layered superconductors on the charge carrier density has been studied in cases of different 
angular momentum of the pairs I, coupling and inter-layer hopping. It has been shown that 
the critical temperature T c is smaller then the 2D critical temperature Tbkt at some values 
of the model parameters, contrary to the XY-model with the parameters J\\ and J± which 
do not depend on carrier density n/, inter-particle coupling V and the temperature of the 
system T. In particular, at small carrier densities T c ^ e F /8, contrary to the dependence 
of Tbkt hi the 2D case. The critical temperature T c is growing with t z , except the case of 
non-zero angular momentum of the pairs at small carrier densities. 

At the same time some questions are remained unresolved. In particular, the behavior 
of the system when the inter-layer coupling t z is not very small has to be studied and the 
doping dependence of the vortex condensation energy should be taken into account. This 
problems are planned to be studied in the future. 
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